Proposed realization of the Dicke-model quantum phase transition 
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The Dicke model consisting of an ensemble of two-state atoms interacting with a single quantized 
mode of the electromagnetic field exhibits a zero-temperature phase transition at a critical value of 
the dipole coupling strength. We propose a scheme based on multilevel atoms and cavity-mediated 
Raman transitions to realise an effective Dicke system operating in the phase transition regime. 
Output light from the cavity carries signatures of the critical behavior which is analyzed for the 
thermodynamic limit where the number of atoms is very large. 
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I. INTRODUCTION 

The interaction of an ensemble of N two-level atoms 
with a single mode of the electromagnetic field is a classic 
problem in quantum optics and continues to provide a 
fascinating avenue of research in a variety of contexts. 
The simplest model of this interaction is provided by the 
Dicke Hamiltonian yj, which takes the form (h = 1) 



H = loo) a + uj Jz + -j= (a f + a) (j + + J_) , (1) 



where u>o is the frequency splitting between the atomic 
levels, ui is the frequency of the field mode, and A is the 
dipole coupling strength. The boson operators {a, a*} 
are annihilation and creation operators for the field, and 
{J±,J Z } are collective atomic operators satisfying angu- 
lar momentum commutation relations, 



2J Z 



J±, Jz 



tJ±- 



(2) 



Contained within the possible solutions to this model are 
a number of significant and topical phenomena, includ- 
ing: 

(i) A zero-temperature phase transition in the ther- 
modynamic limit, N — » oo, occurring at the crit- 
ical coupling strength A c = ^/luloq/2. For larger 
than the critical coupling the system enters a super- 
radiant phase 0, H H |1 H H . 

(ii) An associated change in level statistics, indicat- 
ing a change from "quasi-integrable" to "quantum 
chaotic" behavior 
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(iii) Critical behavior of the atom-field entanglement, 
which diverges at the critical point for N — > oo 

m mum 

It follows that the Dicke model offers a potential setting 
for investigations of quantum critical behavior, quantum 
chaos, and quantum entanglement. 

Practical realization of a system exhibiting the men- 
tioned phenomena presents something of a problem, how- 
ever, in that, in familiar quantum-optical systems, the 
frequencies ui and loq typically exceed the dipole coupling 
strength by many orders of magnitude. This means that 
the counter- rotating terms, a^J + and aJ_ in Eq. l(TJl. 
have very little effect on the dynamics; indeed, they 
are usually neglected in the so-called "rotating-wave ap- 
proximation." Furthermore, dissipation due to atomic 
spontaneous emission and cavity loss is usually unavoid- 
able, significantly altering the pure Hamiltonian evolu- 
tion. Hence, it remains as a challenge to provide a prac- 
tical physical system which might exhibit the interesting 
behavior associated with the idealized Dicke model. 

We propose such a physical system in this paper. We 
suggest a scheme based on interactions in cavity quantum 
electrodynamics (cavity QED) which realizes an effective 
Dicke Hamiltonian Qj| with parameters coq ~ to ~ A that 
are adjustable and can in principle far exceed all dissi- 
pation rates. Our scheme uses cavity-plus-laser medi- 
ated Raman transitions between a pair of stable atomic 
ground states, thereby avoiding spontaneous emission. 
While cavity loss cannot be similarly avoided, it should 
be possible to achieve cavity QED conditions in which 
the dissipation rate from the cavity mode is much less 
than the parameters of the Dicke model. 

In fact, the presence of cavity loss constitutes an im- 
portant and essential aspect of the work presented here: 
output light from the cavity provides a readily measur- 
able signal from which an experimenter can learn, rather 
directly, about the properties of the system. In partic- 
ular, various spectral measurements made on the out- 
put light clearly reveal the critical behavior of the Dicke 
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FIG. 1: Atomic level scheme. Excited states have energies 
hu)j (j = r,s). Such a scheme might be realized, e.g., by alkali 
atoms, with |0) and |1) as different ground-state sublevels. 
Note that |r) and |s) may be the same level, provided that 
the Raman channels remain distinct (which requires uj\ 7^ 0) . 



model as the coupling parameter is changed. 

We begin in Sec. II with a description of the proposed 
scheme for realizing the Dicke model in an optical cavity 
QED system. In Sec. Ill, we briefly discuss a possible 
experimental scenario involving atoms confined within 
a ring cavity and establish parameter values for use in 
the numerical calculations. Our theoretical study of the 
dissipative Dicke model in the thermodynamic limit is 
presented in Sec. IV. It is based upon a linearized anal- 
ysis in the Holstein-Primakoff representation of the col- 
lective atomic spin and the input-output theory of open 
quantum-optical systems. We present results for the cav- 
ity fluorescence spectrum, the probe transmission spec- 
trum, and the spectra of quadrature fluctuations — i.e., 
homodyne, or squeezing spectra. These spectra vividly 



illustrate the changing nature of the system through the 
critical region of the phase transition. We also describe 
a means of computing variance-based measures of atom- 
field entanglement from homodyne spectra of the cavity 
output field. We finish in Sec. V with the conclusion and 
a discussion of possible further investigations. 



II. PROPOSED REALIZATION: BALANCED 
RAMAN CHANNELS 

We consider an ensemble of N atoms coupled simulta- 
neously to the quantized field of an optical cavity mode 
and the classical field of a pair of lasers. All fields 
are co-propagating (in the x direction) TEMoo travel- 
ing waves, with beam waists sufficiently broad compared 
to the atomic ensemble that uniform atom-field coupling 
strengths may be assumed. Each atom has two stable 
ground states, |0) and |1), which are coupled through 
a pair of Raman channels, as shown in Fig. ^ specifi- 
cally, the lasers drive ground-to-excited-state transitions 
|1) <-* \r) and |0) <-> \s) with Rabi frequencies £l r and 
il s , respectively, while the cavity mode mediates the 
|r) «-> |0) and \s) <-> |1) transitions, with dipole coupling 
strengths g r and g s . The detunings from the excited 
states are A r and A s , as shown on the figure. 

With the inclusion of spontaneous emission and cavity 
loss, the master equation for the system density operator, 
p sys , is written as 



Psys 



H S y S , Psys ~t~ ^cavPsys ~T" -^sponPsysi (3) 



where H sys is a sum of Hamiltonians: 



H, 



cav — ^cav^ &j 



(4a) 



for the cavity oscillator, 



N 

#at = Wr l r J') ( r il + U s \Sj) (Sj\ + Wi 

3=1 

+ [(Q r /2)e-^* \r 3 ) {l 3 \e lk ^ + (f> s /2)e~ Mlst \ Sj ) (0j\ e lk ^ + H.c.] } , (4b) 
for the driven atoms (H.c. denotes the Hermitian conjugate), and 

N 

Hint = fa Vi) (0.1 a + g s \ 8j ) (M a) + H.c] , (4c) 
i=i 

I 



for the atom-cavity interaction, where w r , u> S} and u>\ 
are atomic frequencies (see Fig. u)i r and u>i s are the 
laser frequencies, and Xj locates the j-th atom in the 
traveling waves, which have wavenumbers k r , fc s , and k 



(where k r ~ k s ~ k). Cavity loss is included through the 
Lindblad term 

£cavPsys = k ( y 2ap sys a' ! - a' ap sys - Psyg^d) , (5) 



3 



and spontaneous emission through the second Lindblad 

term ^sponPsys- 

From this full master equation a simplified equation 
is derived by neglecting spontaneous emission and adia- 



batically eliminating the atomic excited states. We first 
transform to the interaction picture, introducing the uni- 
tary transformation U(t) = exp(— iH^t), with 



N 

H a = (ljis - w'Jata + {( u l r + w i) l r J') ( r il + ^ l s J') ( s il + ^ l 1 ^ i 1 ^} > ( 6 ) 



where is a frequency close to w\, satisfying 

UJls - W;. r = 2uj[. (7) 

Then assuming large detunings of the fields from the ex- 
cited states, 



|A,J > Or.. 

9r,s: k j "cavi 7i 

where 7 the excited state linewidth and 



(8) 



A r = UJ r - (iOlr + Uj[) , A s = LU S - UJls, (9) 

<5cav = W cav - (UJls -Uj[), (10) 

we make the adiabatic elimination and neglect constant 
energy terms to arrive at the simplified master equation 
for the collective coupling of the ground states |0) and 

|1>, 



f> = -i H,p 
with 

H = loo) a + ujq J z + 5a'aJ ; 



£cav/2, 



(11) 



iJ + + 6' J- 



X, 



a T J+ + a J- 



(12) 



where 



JV 



N 



J+ = 5^11^(0^, J.^lOj)^, (13a) 

3 = 1 3=1 
1 * 

j^ 5 Ed 1 i> <h\-\0i)(0j\) (13b) 



3=1 



are collective atomic operators satisfying commutation 
relations J2J), and the remaining parameters of the model 
are defined by 



LO = 


±N(g?/A r + g 2 JA s )+5 cav , 


(14a) 


UJ Q = 


i(^/Ar-^/A s ) + K-^), 


(14b) 


S = 


g 2 r /A r -g 2 JA s , 


(14c) 


\ r = 


±VNg r n r /A r , 


(14d) 


A, = 


±VNg s Q s /A s . 


(14e) 



With these parameters chosen such that 

g 2 /A r = g 2 /A s , g r fl r /A r = g s Q s /A s , (15) 
H is put into the form of the Dicke Hamiltonian Q , 

H = loo) a + lu q J z + — ^= (a + a 1 ) (j+ + J_) , (16) 

with 



lu = Ng 2 /A r + 5 cav , 
UJ = wi - 
A = ±VNg r fl r /A r . 



(17a) 
(17b) 
(17c) 



Hence, we arrive at a realization of the Dicke model with 
parameters that can be controlled through the laser fre- 
quencies and intensities, and where the characteristic en- 
ergy scales are no longer those of optical photons and 
dipole coupling but those associated with light shifts and 
Raman transition rates. 



III. POTENTIAL EXPERIMENTAL SCHEME 

Before proceeding with our theoretical analysis, we 
pause briefly to outline a possible experimental imple- 
mentation of the proposed model. We imagine the en- 
semble of atoms confined inside a ring cavity where it 
interacts with the quantized cavity mode (field operator 
a) as shown in Fig.|2Ja). The cavity mode copropagates 
with the two laser fields (Rabi frequencies f2 r and fi s ) 
through the ensemble as indicated on the figure by the 
dashed line. Quantized inputs and outputs are assumed 
significant through one cavity mirror only — field opera- 
tors Sin and a out in the figure. 

The atomic excitation scheme might be based on an 
F = 1 <-> F' = 1 transition, as occurs, for example, in 
87 Rb. Such a scheme differs slightly from that of Fig. Q 
and is illustrated in Fig. Hfb). The cavity mode is lin- 
early polarized along an axis perpendicular to an applied 
magnetic field of strength B. The magnetic field splits 
the wif — ±1 sublevels of the F = 1 ground state, al- 
lowing for the excitation of the distinct Raman channels 
shown [13j. 
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Parameter values g r /2ir — 50 kHz, k/2it ~ 20 kHz, 
and N ~ 10 6 appear to be practical [Til Tl5| : thus, 
with the choice f2 r /A r = 0.005, one finds an effective 
coupling strength X/2tt = ^^/Ng r fl r /2T:A r ~ 125 kHz. 
This is significantly larger than the decay rate k, plac- 
ing the system firmly in a regime where the Hamilto- 
nian dynamics can be expected to dominate. Note fur- 
ther that, for these parameters, the spontaneous emission 
rate due to off-resonant excitation of the atomic excited 
states is estimated at |(7/27r)(f2 r /A r ) 2 < 40 Hz, where 
7/27T = 6 MHz has been assumed. Finally, the condition 
lu ~ u>o ~ A can be achieved with appropriate choices of 
the laser and cavity mode frequencies, and ground-state 
level shifts of the order of 2tt-10-15 MHz (> 100A) would 
satisfy the requirement for distinct Raman channels. 

The above set of parameters provides just one exam- 
ple of the possibilities, and a wide variety of parameter 
combinations satisfy the requirements of our model. In 
what follows we concentrate in large part, for numeri- 
cal investigations, on the set of (normalized) parameters 
{uj,ujo,k} = {1,1,0.2}. This choice serves to highlight 
the main physical features of the model proposed. 
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FIG. 2: (a) Ring cavity configuration for implementing the 
proposed realization of the Dicke model. Quantized input 
and output fields are denoted by &i n and a out , respectively, 
(b) Possible atomic excitation scheme based upon an F = 
1 <-> F' = 1 atomic transition and a linearly-polarized cavity 
field a. Note that the magnetic field splittings of the Zeeman 
sublevels are not drawn to scale; the detunings of the optical 
fields from the excited atomic states are much larger than the 
ground-state splittings (i.e., |A r , s | 3> |/i-B|). 



IV. ANALYSIS IN THE THERMODYNAMIC 
LIMIT 

We aim to make a theoretical analysis of the Dicke- 
model quantum phase transition in the thermodynamic 
limit, i.e., for N ^> 1. Our starting point is a semiclas- 
sical analysis of the steady state and its bifurcations, to 
which a linearized treatment of quantum fluctuations is 
added using the Holstein-Primakov representation and 
the input-output theory of open quantum systems. 



A. Semiclassical steady states 

Introducing the c-number variables 

a=(a), /3=(J_), w=(J z ), 



(18) 

where a and (3 are the complex field and atomic polar- 
ization amplitudes, respectively, and w is the (real) pop- 
ulation inversion, we examine the semiclassical equations 
of motion 



a 



i-$=(J3 + /?*), 
VN 



$ = —iuj (3 + 2i—=(a + a*)w, 
v N 



(19a) 
(19b) 
(19c) 



These follow from master equation (|llf> . with Hamilto- 
nian (|16|l and cavity damping @ , by neglecting quantum 
fluctuations and imposing the factorization 



((a- 



Jz 



((a + at))^,), 



((a + at) - ((a + at))<(j_- j + )>. 

The semiclassical equations conserve the magnitude of 
pseudo angular momentum, 

w 2 + \(3\ 2 = N 2 /4. (20) 

We use this conservation law and solve Eqs. I|19a|) (|19cl) 

for the steady states, whence a critical value of the cou- 
pling strength occurs at 



A = A c = -y/ (uj a /uj)(K 2 + U) 2 ). 



For A < A c , there are two steady states, 

a ss = (3 SS = 0, w ss = ±N/2, 



(21) 



(22) 



where the states with negative and positive inversion 
are dynamically stable and unstable, respectively. Both 
states become unstable for A > A c , where the new stable 
steady states are 



X 



Vl-A^/A 4 , 



N 



(3 SS - TyVl-^/A 4 , 



N 



X 2 JX 2 



(23a) 
(23b) 
(23c) 
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These quantities are plotted as a function of the coupling 
strength in Fig. EI Note the bifurcation to states of finite 
amplitude and inversion at A = Af. This is the Dickc- 
model quantum phase transition |20jSS0)I3)H)I3 as 
encountered, without fluctuations, in the thermodynamic 
limit. 




these expressions into the Dicke Hamiltonian, we expand 
the resulting expression under the assumption N 3> 1. 
The goal is to achieve a linearization about the semi- 
classical amplitudes derived above; one must therefore 
distinguish between the so-called "normal" (A < A c ) and 
"superradiant" (A > A c ) phases before the expansion is 
made. 

1. Normal phase (X < X c ) 

The semiclassical amplitudes a ss and /3 SS are zero and 
the expansion is made directly upon the operators as they 
appear in Eqs. 124a(l and l|24h|l . This yields the master 
equation 



(25) 



with the Hamiltonian governing fluctuations (omitting 
constant terms) 

H (1) = wa)a + u Q b r b + X(a 1 + a)(S f + b). (26) 



FIG. 3: Steady-state field amplitude, a ss , polarization ampli- 
tude, /3 SS , and atomic inversion, w ss , plotted as a function of 
the coupling strength A, for uj — uio = 1 and k = 0.2. Only 
stable steady states are shown. 



2. Superradiant phase (X > X c ) 

The semiclassical amplitudes a ss and /3 SS are nonzero 
and the expansion of the Hamiltonian is preceded by 
making coherent displacements of a and b, as both 
bosonic modes are macroscopically excited. Specifically, 
defining 



B. Linearized treament of quantum fluctuations in 
the Holstein-Primakoff representation 

In the thermodynamic limit, N ^> 1, the quantum 
fluctuations are small and may be treated in a linearized 
approach. We follow Emary and Brandes U 0, |n| 
and make use of the Holstein-Primakoff representation of 
angular momentum operators 0, . Collective atomic 
operators J+, J_, and J z are expressed in terms of an- 
nihilation and creation operators, b and w, of a single 
bosonic mode: 



J z = Pb - N/2, 



J- = j\, 



(24a) 
(24b) 



from which, using \b,b*\ = 1, the angular momentum 
commutation relations J2J) are recovered. Substituting 



Jk = X 2 JX 2 < 1, 



we make transformations 

a — > C + CK SS , b — > dt- 



y/(N/2)(l+fiy 



(27) 



(28) 



where a ss and (3 SS are given in Eqs. H23a(l and l)23bjl . and c 
and d describe quantum fluctuations about the semiclas- 
sical steady state. We then proceed with the expansion 
to obtain the master equation 



= -i m 



pi 



(29) 



with the Hamiltonian governing fluctuations (omitting 
constant terms) 



m = ^c+^ (i + fi) dU+ ^ (1 " A)(3 .| A) (d+^) a + Am 



2fi 



8/2(1 + A) 



l + ft 



-(J + c)(di +d), 



(30) 
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and 



c'cp - 



(31) 



3. Eigenvalue analysis 



The quadratic Hamiltonians and dissipative Lindblad 
terms above lead to linear equations of motion for the 
expectation values of c and d. We write 



v = Mv, 
where M is a constant matrix and 



(32) 



(33) 



The eigenvalues of M are plotted as a function of cou- 
pling strength in Fig.0]for U) = loq = 1 and k = 0.2, with 
the four eigenvalues grouped into pairs, one pair associ- 
ated with the "photonic" branch and the other with the 
"atomic" branch; the branches are defined by the A — > 
limit of the corresponding eigenstates (or, in fact, the 
A — * oo limit) 0, E| ■ Note that with the nonzero cavity 
damping, there are two coupling strengths of significance 
in addition to A c ; for u> = luq they are 



A' ~ A r 



A" ~ A c + k 2 /16w , 



(34) 



with A' < A c < A". As A -> [A']_ and A -> [A"] + 
the imaginary parts of the eigenvalues on the photonic 
branch go to zero — respectively, as VA' — A from below 
and \/X — X" from above. They remain zero in the in- 
terval A' < A < A". In correspondence, the real parts of 
the eigenvalues split, with the real part of one eigenvalue 
going to zero at the critical coupling strength A c . 

To complement the figure, in the range < A < k/2, 
the eigenvalues are given by (with lu = luq) 



ep= h = -K/2±iJu$- k 2 /4±A, (35a) 



-k/2 ± i^Juo 2 - k 2 /4=f A, 



with 



A = V^(4A 2 -« 2 ), 



(35b) 



(36) 



where both upper or lower signs are to be taken, while 
for k/2 < A < A c , 



ph 



-n/2±iJuj 2 -K 2 /A- A, 



-k/2 ± i^Juo 2 - k 2 /4 + A. 



Thus we see that e 



P h 



0, e 



P h 



-k, and e± 



(37a) 
(37b) 
-k/2± 

i^2u>Q — k 2 /4 as the critical coupling is approached. 

Above the critical point, similarly simple expressions 
cannot be found. We note, however, that for A > A" the 
photonic branch eigenvalues take on nonzero imaginary 
parts once again, and for large A approach — k ± iujq. 
The atomic branch eigenvalues approach ±iuio/p., with a 
rapidly decreasing real part that scales like // 4 = (A c /A) 8 . 



C. Quantum Langevin equations and input-output 
theory 

The equations of motion of the previous sections con- 
cern the "internal" dynamics of the atom-cavity system. 
To probe this dynamics we consider measurements on the 
light leaving the system through the cavity output mir- 
ror. We make use of the standard input-output theory 
of open quantum-optical systems [T]| [T^, |2(|, |2l| , which 
is nicely formulated in terms of quantum Langevin equa- 
tions for system operators: for A < A c , 



Ka + y2Ka,ia(t), (38a) 
(38b) 



plus the adjoint equations, while for A > A c , 

h = -i c,H (2) - kc+ V2K,a in (t), (39a) 
d,# (2) l, (39b) 



d = 



plus the adjoint equations. The operator fii n (t) describes 
the quantum noise injected at the cavity input (Fig. |2J) 
and satisfies the commutation relation 



a ln (t),aj n (t') =s(t-f) 



(40) 
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FIG. 4: Imaginary parts (upper row) and real parts (lower 
row) of the eigenvalues in the linearized Holstein-Primakoff 
representation as a function of the coupling strength A; for 
ill = u!q = 1 and k = 0.2. Solid (dashed) lines correspond to 
the photonic (atomic) branch. The right-hand column mag- 
nifies the view around the transition at A = A c = 0.5099; note 
the splitting (convergence) at A' ~ 0.5050 and A" ~ 0.5124. 
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In addition, for vacuum or coherent state inputs, one has 
the correlations 



(&in(t),&l(t')) = 6(t-t% 

(al(t),a in {t')) = (a in (i),a in (iO) 



0. 



(41a) 
(41b) 



where (A,B) = (AB) - (A)(B). The cavity output field, 
«out(0' is given in terms of the intracavity and cavity 
input fields as 



a ou t(t) = V2na(t) - a in (t), 



(42) 



from which one calculates the output field correlation 
functions and spectra. 



The quantum Langevin equations are linear operator 
equations. For the purpose of computing spectra, they 
are conveniently solved in frequency space by introducing 
the Fourier transforms 



<9(i/) 



l 



e wt O(t)dt, (43a) 
0\-u) = -= \ e lvt O\t)At, (43b) 

V27T J-oo 



where O denotes any one of the operators a, 6, c, d, or 
a- m . In the resonant case, luq — ui, the solutions are: for 
A < A c , 



J 



a{v) = v2k 
1{ V ) 



V - UJ 



{[k - + Mo)] (v 2 - ljq) - 2iuj X 2 } 5 in (i>) - 2iu Q \ 2 a\ n {-v) 
[re - i(v - wo)} [re - i(y + uo)] {v 2 - us 2 ) + 4w 2 A 2 

[a(v)+at(-vj\ , 



and for A > A c , 



'2k, 



{ [re - i(y + q; )] {v 2 - uj 2 /jl 2 ) - 2iw \ 2 fl} a in (v) - 2iu Q \ 2 jl a\ u {-v) 
[re — i(y — uj )] [re — i(u + wq)] {v 2 — 0J 2 /fi 2 ) + 4w 2 A 2 /i 



AAy/2/(l + A) 
i/-«o(l + A)/(2/i) 



(-!/)]. 



(44a) 
(44b) 

(45a) 
(45b) 



D. Entanglement of the atoms and field 

Quantum fluctuations in the linearized treatment are 
Gaussian, and the solutions to the quantum Langevin 
equations can be used to compute their covariances in the 
steady state. For example, the mean intracavity photon 
number for A < A c is given by 

-i poo />OC 

(o + a)a. = — / / (tf(u)a(v'))dvdis', (46) 

2K J-oo J-oo 

where we need the frequency-space equivalents of the in- 
put correlations l|41a|l and 141b(l . i.e., 

{a ta {y),$ a (i/))=8{y- V r), (47a) 
(4»,a in (j/)> - (ainM.ain (!/')) =0. (47b) 

The computed output photon flux from the cavity, 
2n(a) a) ss , is plotted for several different values of re in 
Fig. illustrating a "smoothing-out" of the transition 
with increasing cavity linewidth. The mean excitation of 
the atomic mode, (b j< b) ss , shows similar behavior. 

Of particular interest is the behavior of the bipartite 
quantum entanglement in the vicinity of the critical point 
[3|llllllil|2l|2l|2l|2l|23. The cavity and atomic 



modes are natural choices for the entangled subsystems, 
and given that their fluctuations are described by a Gaus- 
sian continuous variable state, the criterion for insepara- 
bility can be formulated in terms of the variances of ap- 
propriate subsystem operators. In particular, we define 
the quadrature operators 

K = \{ae- ie +a)e ie ), (48a) 

X* = i (be-** + S+e 4 *) , (48b) 

with adjustable phases 9 and <j>, and introduce the EPR 
(Einstcin-Podolsky- Rosen) operators 

u = X 8 a +X*, v = X e a +*/ 2 -X£ +7r/2 . (49) 

Then a sufficient condition for the inseparability of the 
state below the critical point is (any 9 and 4>) 

((An) 2 ) + ((At,) 2 ) < 1. (50) 

Alternatively, a stronger condition may be given in the 
modified form |2^| 

{{Au) 2 ){{Avf) <\. (51) 



8 



4 

3.5 

X 3 

C 2.5 
O 

1 2 
ft 

3 1 5 

B* 

3 

O 1 

0.5 



0.3 

A 



FIG. 5: Output photon flux as a function of coupling strength; 
for lu — ujo = 1 and k = 0.1 (solid), 0.2 (dashed), and 0.5 
(dot-dashed). 



the probability of detecting a photon of frequency v at 
the cavity output, (ii) the probe transmission spectrum, 
the transmitted intensity as a function of frequency of a 
(weak) probe field applied at the cavity input, and (iii) 




Above the critical point, similar definitions and condi- 
tions based on operators c and d hold. Here the EPR 
variance is associated with a quantum state "localized" 
about one of the two possible semiclassical steady states 
(|23a|l - H23c|l : within our linearized treatment transitions 
between these states are ignored. 

The sum of EPR operator variances — inequality (|50J) — 
is plotted as a function of the coupling strength in Fig. [15] 
It approaches a cusp-like minimum at the critical cou- 
pling strength; thus, the entanglement here is maximum. 
The variance product — inequality l|51l) — exhibits similar 
behavior. These variances are measurable quantities. 
They offer a means of tracking entanglement across the 
phase transition. In fact, as we show in Sec. IIV E 41 
variance-based entanglement measures can, under appro- 
priate conditions, be deduced from measurements on the 
cavity output field alone. 



E. Spectra of the cavity output field 

Cavity output field spectra can be computed from the 
solutions to the Langevin equations (|44a(l and (|45a|l . the 
correlations (|47a|l and il47b(l . and the input-output rela- 
tions 



Oout(^) = V2na(is) - a la (v), 



(52) 



A < A c , and 
a ou t 0) = 



'2k, 



c(v) + \f2~Ka ss 8(v) -a in (is), (53) 



A > A c . We consider three standard spectra: (i) the flu- 
orescence (or power) spectrum, which is proportional to 



FIG. 6: Sum of EPR operator variances as a function of cou- 
pling strength; for u = a>o = 1, K = 0.2, 6 = tan - (k/uj), and 
<jf> = 0. The choice of 9 minimizes the sum in the vicinity of 
the critical coupling strength, A = A c — 0.51. 



homodyne spectra, which measure the quantum noise 
variances of output field quadrature amplitudes. 

1. Fluorescence spectrum 

The fluorescence spectrum consists of a coherent part, 
representing the mean excitation of the intracavity field, 
the semiclassical solution a ss , and an incoherent part 
which accounts for the quantum fluctuations. The lat- 
ter is defined by 



(5out M> a out (>')) = S{v)8(v - v'). 



(54) 



It can also be expressed in terms of the steady state au- 
tocorrelation function of the intracavity field, with 



A < A c , and 



e-^(at( T ),a(0)) ss dr, (55) 



<ct(r),c(0)) ss dr, (56) 



A > A c . Making use of solutions 144a|) and (|45a|l for a(v) 
and c(i^), and the input correlations 147a|) and l|47b|) . one 
finds 
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4kw A 2 m (1 ' 2) 



k - 1(1/ - wo)][« - i(f + wo)] ^ 2 - w 2 / (/x^ 2 )) 2 + 4w 2 A 2 /x( 1 > 2 ) 



(57) 



with the definitions 

= 1, m (2) = A c 2 /A 2 . (58) 




FIG. 7: Incoherent part of the cavity fluorescence spectrum 
S(u) for various values of coupling strength A; for w = u>o = 1 
and k = 0.2 (A c ~ 0.51). 

Sample spectra S^ 1,2 ^(u) are plotted in Fig.[7| The po- 
sitions and widths of the spectral peaks are determined 
by the eigenvalues of the linearized dynamics discussed in 
Sec. IIV B 31 Thus, below the critical point the spectrum 
shows central and outer doublets associated with the pho- 
tonic and atomic branch eigenvalues, respectively. The 
peaks of the photonic branch doublet merge as A — > A c , 
forming a single narrow peak at v = 0; within the lin- 
earized treatment the intensity under this peak diverges 



at A = A c . Above the critical point a pair of doublets 
appears again. Far above the critical point the pho- 
tonic branch peaks approach detunings determined by 
the cavity mode resonance frequency, v ~ ±lu = ±1, 
and linewidths (FWHM) 2k = 0.4. The atomic branch 
peaks move linearly apart, following the increasing Rabi 
frequency in the presence of the increasing mean intra- 
cavity field; they also become increasingly sharp. 

Note that the symmetry of the spectra is ensured 
by energy conservation and the fact that, due to the 
symmetrical nature of the atom-cavity coupling, photon 
emissions from the cavity can be associated with transi- 
tions to either lower or higher internal energy states of 
the atom-cavity system. 

2. Probe transmission spectrum 

One may also examine the system by driving the cavity 
mode with a (weak) laser field and measuring the inten- 
sity of the coherent transmission as a function of laser 
frequency. Such a measurement provides a rather direct 
probe of the energy level structure of the atom-cavity 
system; only when the laser frequency matches a system 
resonance would substantial transmission be expected. 

Analytically, we treat the measurement by adding a 
driving term, £ p e~ Wpt , to the equations of motion for a 
and c, where £ p and v v are the probe field amplitude and 
frequency. Solving the equations of motion in frequency 
space as before, the coherent amplitude in transmission 
follows straightforwardly from the coefficient of 8(v — v p ) 
in the solution for (a ou t(f))- The transmitted probe in- 
tensity is thus found to be 



[k - i{v v + loq)} 






[k - i(v p - ojq)] [k - i{is p + wo)] 




+ 4u; 2 A 2 £< 1 ' 2 ) 



(59) 



with fl^ and p,^ defined by Eq. flBl^l : the normalization 
is such that the spectrum is a Lorentzian of width 2k and 
unit height when A is set to zero. 

A series of probe transmission spectra are plotted in 
Fig. [SI where we choose values of coupling strength to 
correspond to Fig. {7\ The spectra contain two principal 



peaks, one associated with the photonic and one with 
the atomic branch. Their behavior as a function of A 
replicates the behavior displayed by the fluorescence. 
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FIG. 8: Probe transmission spectrum T(v p ) for various values 
of the coupling strength A; for w = u>o = 1 and ft = 0.2 
(Ac ~ 0.51). 
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3. Homodyne spectra 

Homodyne spectra measure the fluctuation variances 
in frequency space of the output field quadrature am- 
plitudes. Quadrature operators are defined in time and 
frequency space, respectively, as 



X ou t,e(v) 



-w 



J8 



2 

1 r 

2 



a ut {v)e 10 + 5ol 



(60a) 
(60b) 



where 9 is the quadrature phase. The (normally-ordered) 
ho mody ne spectrum, S out g{v), is defined by the variance 



out,e(v),X ou t,e(v') ■) = S oll t,o(v)6(v + i/), (61) 

which we compute from the input-output relation l|42(l 
and solutions l)44all and (|45a|l for the intracavity fields. 
Note that with the choice of normal ordering the vac- 
uum noise level corresponds to S ou t,e(v) = 0, while per- 
fect quantum noise reduction corresponds to S ou t e {v) = 
-1/4. 

Numerical results for 9 = and 9 — tt/2 are pre- 
sented in Fig. As the coupling strength approaches 
A c , the phase transition is signaled by a divergence of 
the quadrature amplitude flutuations at v — 0, simi- 
lar to the behavior of the cavity fluorescence (Fig. [7J). 
Nonetheless, there is an optimal 9 at each A, for which 
near-perfect noise reduction occurs in the (9 + tt/2)- 
quadrature at v — 0. Figure HUI plots the optimal phase 
and corresponding minimum quadrature variance across 
the threshold region. As A — > A c , the optimal phase ap- 
proaches 9 m m = t&aT 1 {k / lo) + n/2. 

Well above the critical point the noise level returns to 
the vacuum noise level at all frequencies, except close to 
the atomic branch resonances at v ~ ±u>o/jl. Here sig- 
nificant squeezing below the vacuum noise level is found 



FIG. 9: Quadrature noise spectra S out> e(v) with 8 = (solid) 
and 8 = n/2 (dashed); for u) = u>o = 1> K = 0.2 (A c ~ 0.51), 
and A = 0.4, 0.49, 0.505, 0.52, 0.6 (top to bottom). 



for the 9 — quadrature amplitude, with correspond- 
ing amplification of the fluctuations at 9 = tt/2. In fact, 
substantial squeezing of the atomic branch resonances oc- 
curs also for A < A c , as seen from Fig. |5J Although the 
bandwidth of this squeezing becomes increasingly narrow 
as A increases, the noise reduction on resonance actually 
approaches 100%. 




0.3 0.35 0.4 0.45 0.5 0.55 0.6 0.65 0.7 




FIG. 10: Optimal squeezing at v — as a function of A (top) 
and the quadrature phase angle at which the optimum occurs 
(bottom); for ui = uo = 1 and n = 0.2 (A c — 0.51). 
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4- Output field squeezing and atom-field entanglement 

For the parameter regime we have considered, the spec- 
tra presented exhibit distinct features that can be iden- 
tified with either the "photonic" or "atomic" branches. 
The lower frequency peaks are associated with the pho- 
tonic branch and the higher frequency peaks with the 
atomic branch. The corresponding photonic and atomic 
modes are formalized by a diagonalization of Hamilto- 
nians l|26|l and (|30l) via Bogoliubov transformations, as 



shown in [9j and outlined in Appendix A. If these modes 
are well separated in frequency, and k is sufficiently small, 
then one can also associate with them what are essen- 
tially independent and uncorrelated output fields, a a ut (i) 
and ciout (^); hence, we can relate their quadrature vari- 
ances to the quadrature variances of linear combinations 
of the "bare" internal atomic and cavity modes. In par- 
ticular, in the normal phase, we find that the EPR vari- 
ance of Eq. (|50|l (with <f> = 6) is approximately given by 
(Appendix A) 



- (( : ^out,0: '■) + (■ ^out,0+7r/2' + ^ ( 62 ) 

where the output field quadrature variances are calculated from integrals of the (normally-ordered) homodyne spec- 
trum over appropriate frequency ranges, i.e., 

( : ^out.e+V2'^out,e+V2 : ) = 5= / Sout,e+*/2(v)dv, (63a) 
(:A a uM ,A a uM :> = i- f S ont , e (u)du. (63b) 

J { ^at } 

If one then considers the homodyne spectra plotted for A = 0.4 and 0.49 in Fig. |5| qualitatively, these expressions 
allow entanglement to be inferred from the fact that S'out,o(i / ) exhibits squeezing — i.e., is negative — on the atomic 
branch while S ont ^i 2 (v) exhibits squeezing on the photonic branch. Given the well-defined peaks and dips in the 
homodyne spectra around 8 = 0, ir/2, we estimate <|62[) by 

V est = -^-< / S , out,e(^)di'+ / 5'out,e+7r/2(^)d^ > + 1. (64) 

K27T {■> Saut, 9 O)<0 Js outi8 + „ /2 H<o J 

I 



This quantity is plotted as a function of A in Fig. ^2 
For A < A c , it shows rather good agreement with the 
EPR variance plotted in Fig. [5J the agreement improves 
for decreasing values of the decay rate k. Above thresh- 
old, on the other hand, V cst can only be regarded as a 
good measure of entanglement when A is quite close to 
A c . In the superradiant phase, the relationships between 
output field and internal mode operators are more com- 
plicated [compare Eqs. (|A6a|) - (|A7bjl and (|A13aj) - l|A17bj . 
but entanglement measures based on output field quadra- 
ture variances can still be derived. The measures depend 
explicitly on A, however, and cannot be directly related 
to the EPR variance of Eq. H5U|I , as was possible for the 
normal phase. Nevertheless, they do display a drop-off in 
the degree of entanglement with increasing A, consistent 
with that shown in Fig. 



V. CONCLUSION 

In this paper we have proposed a scheme for the real- 
ization of an elementary atom-light interaction Hamilto- 
nian - the so-called Dicke Model - which should enable 



the observation and detailed study of a quantum phase 
transition involving a collective atomic pseudo-spin and a 
single quantized mode of the electromagnetic field. While 
the optical cavity-QED system considered is necessarily 
dissipative, due to cavity loss, the dissipation is a positive 
feature providing a window through which one can moni- 
tor the system using standard quantum-optical measure- 
ment techniques. As we have demonstrated, fluorescence, 
probe transmission, and squeezing spectra all provide de- 
tailed information on the varying energy level structure 
of the Dicke Hamiltonian and exhibit striking behavior 
in the vicinity of the critical point. 

We have focussed exclusively on the thermodynamic 
limit, with the number of atoms taken to infinity, where 
fluctuations may be treated using a bosonic approxima- 
tion for the collective atomic spin and linearization about 
the semiclassical steady state. Finite-size systems are a 
natural consideration, both theoretically and experimen- 
tally, and are of interest for examining scaling properties 
and deviations from the linearized model. Indeed, in a 
regime of strong-coupling cavity QED (see, for example, 
[30l I3TL 13^ 1 it could be possible to realize the critical 
regime of the Dicke Model with just a few atoms. In 
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FIG. 11: Estimate of the EPR variance, Vest, as a function of 
coupling strength; for u = uio = 1, K = 0.2, and quadarture 
phase angle 6 = tan - (n/ui) — 0.2. The choice of 6 minimizes 
Vast close to the threshold at A c ~ 0.51. 



(possibly including other cavity modes) one could en- 
visage making independent measurements on the atomic 
ensemble to complement those made on the cavity field, 
enabling the explicit determination of correlations and 
entanglement measures such as the EPR variance. Sep- 
arately addressable atomic sub-ensembles coupled to the 
same quantized cavity mode would also allow the mea- 
surement of entanglement between different "blocks" of 
spins [27^ 
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such a case, issues of quantum measurement (e.g., mea- 
surement backaction) arise, providing a further interest- 
ing avenue of investigation. 

Finite-size systems and a full treatment of the Dickc 
model without linearization are also of importance for 
studying the role of quantum entanglement in the vicinity 
of the phase transition. In the present paper we touched 
only briefly on this subject, demonstrating that variance- 
based measures of atom-field entanglement can in princi- 
ple be determined from homodyne spectra of the cavity 
output field, thus enabling entanglement to be "moni- 
tored" . The proposed system clearly offers further excit- 
ing prospects for the study of entanglement in a quantum 
critical system. For example, with additional light fields 



APPENDIX A: NORMAL MODES AND 
ENTANGLEMENT CRITERIA 

1. Normal phase 

The normal-phase Hamiltonian (|26[) can be diago- 
nalised in the form (omitting constant terms) 

=u$(\)AlA + u>$(\)&B, (Al) 

where w^ (A) and u^P(A) are the normal mode frequen- 
cies, with respective normal mode operators [t| 



A = \ (iJu^g) - wo)@ - St) + (0$ + cj )(a - b) 



-1/2 



-1/2 



where u) — loq has been assumed. The inverse relationship for the cavity mode operator a is 



1 



& = 2 i l 2uJaUJ ph 



(i) 



-1/2 



( Wo - W W)it + ( Wo +a,W)i 



ph 



-1/2 

2cJ Wat 



(A2a) 
(A2b) 

(A3) 



If the normal modes are well separated in frequency with 
linewidths much smaller than their separation, within the 
bandwidth of the photonic mode the cavity mode contri- 
bution to the output field [Eq. (35)1 ] may be written as 



1 

2 V 



-1/2 r 



(wo 



4^ + (^ + ^)A 



(A4) 
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and within the bandwidth of the atomic mode as 



-1/2 



(A5) 

Using these approximations, the input-output relation 
(|42p. and Eqs. l|A2afl and (|A2 b|) . one may derive approx- 
imate expressions for the output field quadrature opera- 
tors in the specified frequency regions in terms of "bare" 
cavity and atomic mode operators: 



ph 

out,e 



^(X e a +X e b )/2-X?} e 



(A6a) 
(A6b) 



It follows that, in the normal phase, the normally-ordered 
output field variances can be directly related to the in- 
ternal mode EPR variances [l£j : 



:X 



P h 



out,6>' out,S 



■ ^-out,0 



,xf 



out,e ■ 



:X. 



Y t) Vb 



v yfl 



xl. 



(A7a) 
(A7b) 



where a vacuum field input has been assumed. Then, 
adopting the entanglement criterion from 28], entangle- 



ment between the cavity and atomic modes can be in- 
ferred whenever the inequality 



' ^out,e+7r/2' ^out,0+7r/2 : ) < 



is satisfied. 



(A8) 



2. Superradiant phase 

The superradiant-phase Hamiltonian (|3(J|) can be diag- 
onalised in similar fashion in the form (omitting constant 
terms) 



(A9) 



(2) (2) 

where Wp h and ui at are the above threshold normal mode 
frequencies (for linearization around cither of the above 
threshold steady states) , and the respective normal mode 
operators are given by the somewhat more complicated 
expressions 



J 




D = 



P h 



U) )£ 



sm7 



(2) 



u )c? + (ul t > + w )c 



(2) 



cos 7 



(2) 



(2) 



(wg;-cDo)(? + (wS+Oo)d" 



d / 2) 



ph 



ujo)^ + (w^ 2) + w )d 



(AlOa) 
(AlOb) 



r 



with 



tan(2 7 (2) ) = 2/t 2 (l-/i 2 )- 



(All) 



pic a relationship between output field and internal mode 
quadrature variances to be written down. Nevertheless, 
following the same arguments as before, one can write 



loq = uo(l + ft 1 )/2 



(A12) 



where, once again, the resonance condition to = luq has 
been assumed. These expressions do not allow for as sim- 



where 



^out,6> 

Xat 
out,e 



'2n XI 



X\ 



P h 



<2kY, 



Cll 



X 



in, 9i 



(A13a) 
(A13b) 



J 



; os 2 ( 7 ( 2 ))X c e -cos( 7 ( 2 ')sin( 7 ( 2 )) 
^ = sin 2 ( 7 ( 2 ))^ + cos( 7 ( 2 ))sin( 7 ( 2 )) 



ILVq - 6I=-it/2 
77 d, 



a r°Q v0=O , • n U ° v e =K/2 

costM — X d +sm£M — X, 1 
V w V w 



(Al4a) 



(A14b) 



For these more complicated linear superpositions of the internal mode operators it is still possible to derive insepara- 
bility criteria based on their variances. In particular, following |29j one can show that a sufficient condition for the 
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inseparability of the system state is given by 

V x = 



(A^ d ) 



cos 2 ( 7 ( 2 ))sin 2 ( 7 ( 2 )) 



< 1, 



(A15) 



or in stronger form 



_((AxL + ' /a ) 2 )((Ay^ 

±cos 4 ( 7 ( 2 ))sin 4 ( 7 ( 2 )) 



< 1. 



(A16) 



The required variances can be deduced from the (normally-ordered) photonic and atomic output field quadrature 
variances by inverting the relations 



c<»8-(7«) + cos 2 ( 7 ( 2 )) sin 2 ( 7 ( 2 )) ( ^ sin 2 9 + ^ cos 2 i 

LU LUq 



and 



2kU(AY c 



e \2 

cd) 



r 4 ( 7 ( 2 )) +cos 2 ( 7 ( 2 ))sin 2 ( 7 ( 2 )) ^ cos' + — sin' 

\U>0 ^0 



(Al7a) 



(A17b) 



Numerical examples of V% and V2 versus A are shown in Fig. 1121 The output field quadrature variances used were 
computed via numerical integration of the homodyne spectra from Section IV. E. 4. The computed V% and V2 display 
a decay in the degree of entanglement with increasing A, consistent with the internal mode EPR variance of Fig. 




0.75 



FIG. 12: Entanglement measures Vi (dashed) and V2 (solid) versus coupling strength; for ui = u>o = 1, k = 0.2, and — 0. 
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